Characterization of coherent impurity effects in solid state qubits 
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We propose a characterization of the effects of bistable coherent impurities in solid state qubits. 
We introduce an effective impurity description in terms of a tunable spin-boson environment and 
solve the dynamics for the qubit coherences. The dominant rate characterizing the asymptotic time 
limit is identified and signatures of non-Gaussian behavior of the quantum impurity at intermediate 
times are pointed out. An alternative perspective considering the qubit as a measurement device 
for the spin-boson impurity is proposed. 
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Coherent nanodevices are inevitably exposed to fluc- 
tuations due to the solid-state environment. Well stud- 
ied examples are charged impurities and stray flux tubes 
which are sources of telegraphic noise in a wide class of 
metallic devices. Large amplitude low-frequency (mostly 
1//) noise, ubiquitous in amorphous materials is 
also routinely measured in single-electron-tunneling de- 
vices [2j]. Noise sources are sets of impurities located in 
the oxides and in the substrate, each producing a bistable 
stray polarization. Telegraphic noise has also been 
observed in semiconductor and superconductor based 
nanocircuits Q . The possible presence of impurities en- 
tangled with the device has been suggested in [4j . Recent 
experiments on Josephson qubits indicated that charged 
impurities may also be responsible for noise [f| exhibiting 
an ohmic power spectrum at GHz-frequencies. Different 
theoretical models have been proposed aiming to a uni- 
fied description of broadband noise sources. They share 
the common idea that the variety of observed features 
are due to the dynamics of ensembles of bistable impuri- 
ties [1, @, 0, @| ■ In particular in Ref. [§] it has been pro- 
posed that a noise power spectrum compatible with the 
observed relaxation of charge- Josephson qubits H can be 
obtained if sets of coherent impurities are considered. 

Solid-state noise also determines dephasing. This issue 
has attracted a great deal of interest in recent years since 
it has been recognized as a severe hindrance for the imple- 
mentation of quantum hardware in the solid state. The 
effect of slow noise due to ensembles of thermal 0, Hoj 
and non-thermal Q fluctuators has been addressed. Slow 
noise explains the non-exponential suppression of coher- 
ent oscillations observed when repeated measurements 
are performed 
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Il2j . In addition fluctuations active 
during time evolution represent an unavoidable limitation 
even when a single-shot measurement scheme or dynam- 
ical decoupling protocols [T3| are available. Note that at 
experimental temperatures (~ 10 mK) quantum impuri- 
ties may have a significant influence. 



In this Communication we investigate qubit dephas- 
ing during time evolution due to coupling to a coher- 
ent impurity. The full qubit dynamics is solved in the 
regime where qubit relaxation processes are absent. We 
show how the coherent and non-linear dynamics of the 
impurity is reflected in the qubit behavior. We identify 
regimes characterized by strong qubit - impurity back- 
action. Specifically, we discuss dependence on the impu- 
rity preparation and beating phenomena. An alternative 
interpretation with the qubit acting as a measurement 
device for the impurity is presented at the end of this 
Communication. 



Modeh 
tem, "Hi 
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We model the impurity as a two-state sys- 
coupled to the qubit (a) 
This anisotropic coupling 
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(71 = 1). 



has been discussed for charge qubits, where it models 
the electrostatic interaction [H, 0, [l(| • In this case the 
two physical states (r z — * ±1) correspond to a bistable 
stray polarization of the qubit. They are viewed as the 
ground states of a double- well deformation potential, the 
impurity oscillating coherently between them with fre- 
quency £li — \J e 1 + A 2 . Dissipative transitions between 
the minima come from the interaction with a bosonic 
bath 



m {Kb = E c 



via Hib 



The 



operator X = J2 a ^a(a a +a' a ) is a collective displacement 
with ohmic power spectrum S(u>) = 2tt Klo coth ^ with 
a high-energy cutoff at uj c (k B — 1) (l4| . This spin-boson 
environment (SBE) may induce a variety of qubit dynam- 
ical behaviors, since its deg ree of coherence depends on K 
and on temperature T LJ] . For instance for weak damp- 
ing, K -C 1 a crossover occurs between a low "impurity 
temperature" , T -C VLi regime, where the impurity per- 
forms damped oscillations, to the regime of incoherent 
dynamics if T > tti (white noise S(u>) » AirKT) 

We assume that the qubit Hamiltonian conserves o~ z , 
therefore the impurity induces pure dephasing [3] with 
no relaxation of the qubit [l6| . This regime is very inter- 
esting since energy exchange processes do not blur deco- 
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FIG. 1: (a) Impurity Bloch sphere. An isolated impurity 
Hi defines the mixing angle 6 = arctanA/e, Ti± define 
e± = arctanA/(e ± v). (b) Impurity bands ±^/E 2 + A 2 : 
impurity energy splittings depend on the qubit state, fl± — 
\/ (e ± v) 2 + A 2 . Eigenstates of TLo, are { | i } }, i = a, b, c, d. 
Conservation of a z allows only intra-doublet processes, a <-> b, 
c <-> d. 

herence of the qubit, which is then maximally sensitive to 
the SBE dynamics. Pure dephasing due to Fano impuri- 
ties was addressed in [9|, recently the asymptotic dynam- 
ics has been studied [17] . This model corresponds to a 
over-damped impurity (SBE at K = i), here we consider 
K <C 1 where the impurity may behave coherently. 

Method and analytic results. — For pure dephasing the 
qubit Hamiltonian can be gauged away by a proper ro- 
tation. In this picture we consider the reduced density 
matrix p{t) describing the entangled qubit-impurity sys- 
tem. For K <C 1 the interaction with the bosonic bath is 
studied by the Born-Markov master equation (ME) [l8| 

d tP (t) = -i[H , P (t)] - Jltt'[\S{t>) [T Z , [T z {t'), P {t)\] 

+ UV)[t„ [r»(tV(*)] + ]}, (1) 

where 7io = Hqi + Hi is the undamped Hamiltonian. 
Here, the transform S(t) of the power spectrum and the 
bath susceptibility = -i{[X(t),X(0)]+)Q(t) enter 
the damping term. We introduce the conditional Hamil- 
tonians of the impurity TL± = — ^ (e ± v) t z — | A t x , (see 
Fig. [T]) and the eigenvectors of Ho, I \i ) \, which are fac- 
torized in eigenstates of a z and of Ti± [15[ . The qubit dy- 
namics at pure dephasing is described by the coherences 
(o~±(t)) = Tr (a x ± ia y ) ® l r ], and in particular 

(<r_(i)) = [/?ac(i) + /»bd(i)] cos0+ [/?<,<*(*) — /3& c (t)] sia0, 

where = |(#- — 9+) is a combination of the mixing an- 
gles of H.± (Fig. [T]). Since a z is conserved, the damping 
tensor presents only four non vanishing 4x4 diagonal 
blocks. We focus on the block acting on the terms enter- 
ing (cr_(i)). Performing a partial secular approximation 
within this block, we get two sets of decoupled equations 
for p aci ptd and p a d, Pbc- We quote here the first set 

pu(t)) { r 21 -i5-t 2 ) {p M (t))' 

where S = j(^+ — Fig. [1] The rates Ti, describ- 

ing dissipative transitions and pure dephasing processes 



between the 4-states and the bosonic bath, read 

ri, 2 = a+r T (^+) + a 2 _r T (fi_) + Vs s(o) , 
r 12i2 i = a + a_[r ± (fi + ) + r ± (fi_)], (3) 

a± = ^=sin0±; rj s = \ sin 2 6> sin 2 0, 

where 6 = \(9+ + 0_). Here T±(u>) = 2vrA^[coth(^) ± 
1], are the impurity emission (+) and absorption (— ) 
rates of energy u>. The elements p a d, Pbc satisfy similar 
equations with S replaced by fl = j(^+ + fi-) and rates 

T 3 ,4 =alr T (n + )+a 2 _T ± (n_)+r ]c S{0), 

r 3 4, 43 = a+a- [r± (n+ ) + r T (n_ )] , 

7 h = I cos 2 9 cos 2 4> . 

Diagonalization of Eq. and of the corresponding set 
f° r Pad, Pbc yields the eigenvalues 

Ai. 2 - -^f 2 - ± \y/{2iS + r 2 - ro 2 + 4r 12 r 21 , 
a 3 , 4 = - £a ^ ± ± v / ( 2 ^ + r 4-r 3 ) 2 + 4r 34 r 43 . 

The explicit form of (<r_ (t)) depends on the initial con- 
ditions for pit). Because of the high accuracy of prepa- 
ration presently achieved in solid state implementations, 
factorized qubit-impurity states p(0) = Pcr(O) ® /O T (0), 
represent a realistic scenario. The impurity initial state is 
instead out of the experimentalist control, thus we choose 
p T (0) = | (1 T + p z t z ), p z being the initial average of t z . 
The impurity starts from a totally unpolarized state for 
p z = 0, from a pure state if p z — ±1. This class of initial 
states guarantees the positivity of the dynamical process 
ensuing from Eq.([T]). With this choice we find 

( ? _(t)) - (MO)} E^ 4 , (5) 
A ii = 2(xf^) i C0S ^ t A i ~ X2± ( r i2 + r 2 i)]T 
Pz cos{9 + 9)[-2iS -r 2 + r 1 + r 12 - r 21 ]} , 

(6) 

^3,4= 2 j0 x ^ {sin0 [A3 - A 4 T (r 34 + T 43 )]± 

Pz sin(9 + 9) [-2*n - r 4 + r 3 - r 34 + r 43 ] } . 

Eqs. (J4j)- ((6]) are the main result of this Communication. 
They cover the parameters regime where S(£l±) <C Cl±. 
Single-phonon processes dominate at low T, whereas 
multiphonon-exchanges are paramount at higher T where 
the white noise results of |15[ are recovered. Reliability of 
ME is confirmed by a real-time path-integral calculation. 

Discussion of the results. — We focus our analysis on 
the low-temperature regime T <C f2_. Here effects of 
the dissipative processes internal to the SBE on the 
qubit behavior are clearly identifiable. In this limit en- 
ergy absorption processes are exponentially suppressed 
(r_(f2±) w 0) and the eigenvalues take the forms 

Ai =iS- Vs S(0), 

A2 = _ + 7°+ + 7r- + 7°- _ ; (?) 

A 3i4 = ± ifi - 7rT ± - TfcS(O) , 
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FIG. 2: The four rates TZe[Xi] from Eqs.©. In (a) T = 0. 
Inset: imaginary parts (independent on temperature for T < 
In (b) T = 0.5A. Parameters are K = 0.1, e = 3A. 



where intra-doublet relaxation rates (see Fig. [T]) 

7r± = ^sm 2 (0 ± )5(O ± ) = I(A) 2 5(fi±) (8) 

have been introduced (7° ± value at T — 0). Note that 
pure dephasing processes oc 5(0), are not simple sum of 
intra-doublet dephasing terms, 7^± = i cos 2 (6* ±) 5(0). 

In the following we present a selection of illustrative 
behaviors for e > A. In this regime the two conditional 
Hamiltonians TL± may differ significantly and enforce pe- 
culiar impurity dynamical behaviors. For example beat- 
ings when 8 approaches f2, i.e. around e — v which iden- 
tifies a sort of "resonance regime" for our problem. 

We first characterize the asymptotic qubit dynamics, 
by the T and v dependence of the eigenvalues. At zero 
temperature the pure dephasing contributions fade away, 
and one rate, TZe [Ai], vanishes, as expected. Only emis- 
sion processes contribute to the residual rates, and they 
directly sound out intra-doublet relaxation rates y®±. 
Their behaviors reflect the sensitivity of Tt±, to noise 
acting along r z . While 7"+ decreases with increasing v, 
7°_ takes a maximum at the resonance point (see Eq.flSJ), 
the "transverse" (0_ = tt/2) noise condition for 7Y_. This 
implies a non-monotonous dependence oiTZe [A2,a] on the 
coupling v, Fig. [2ja) . The imaginary parts of Ai,2 and 
A3_4 interchange characters at resonance (Fig. inset) 
leading to possible hybridization (see below). Increasing 
T leading correction to the rates come from pure dephas- 
ing terms 5(0). As a difference with T = 0, all the rates 
are finite and cross around resonance, Fig. [^b). 

These features are crucial for the asymptotic dynamics 
of (er_ (i)) , which does not depend on the impurity prepa- 
ration. We then expect at T = 0, undamped oscillations 
with 8, while at finite T, damped oscillations driven by 
one or two complex eigenvalues. For example in the case 
of Fig. [2](b) the dominant rate is TZe [Ai] for v < e and 
TZe [A4] for v > e. It is a non-monotonous function of v 
and a cusp signals crossing of eigenvalues (a similar effect 
may explain non-monotonic behavior of 171]). 

At intermediate times, all eigenvalues may be relevant, 
depending on the weights Ai in Eq.(j6|). To substantiate 
this point, in Fig. [3] we show Ai corresponding to the 
eigenvalues in Fig. [^a) for different preparations. Re- 
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FIG. 3: Weights \Ai\ of e A *' from Eq.© as a function of v/e. 
(a) Dominant weights in the small v region: |Ai| (full sym- 
bols) and I j4.2 | (open symbols) for p z = (squares), p z — — 1 
(circles) and p z = 1 (diamonds). Effect of impurity prepa- 
rations: p z = 1 (b), p z — —1 (c), and p z = (unpolarized 
state) (d). I Ai| (blue), |A 2 | (blue dashed), |A 3 | (red), \A 4 \ 
(red dashed). Parameters: T = 0, e = 3A, K = 0.1. 



markably, the weights very weakly depend on T (not 
shown), then the following picture generally holds for 
T < For extreme weak coupling, v — > 0, \A\\ w 1 
(Fig. EJa)) implying universal dynamics independent on 
the initial conditions. The dominant eigenvalue is Ai 
with S — > and (<j-(t)) decays exponentially with the 

Golden rule rate Tor — \ jt+s 5 ' sm4 ^' ^ n ^ ms re gi me 
the impurity acts as a Gaussian reservoir and may be 
described with linear response theory in the coupling 
v. Away from this tiny region non-Gaussian effects oc- 
cur and different impurity preparations result in different 
time behaviors, giving separate information on the vari- 
ous eigenvalues. Far from resonance, a single frequency 
shows up in (<r_(i)) independently on p z (5 if v < e, O 
if v > e). Damping of the oscillations depends on the 
initial condition, Fig. [3] (b) - (d). For instance, at finite 
v < e, the decay occurs with IZe [X±] if p z = 1 and with 
IZe [A2] if p z = — 1) both rates are present for unpolarized 
initial state. This behavior is stable against temperature 
variations. Beatings and T - dependence are instead char- 
acteristic of the resonant regime. At v = e, at least two 
amplitudes are equal, |Aj.| w \A^\ (p z = 1) or \ A 2 \ ~ |A 3 | 
(p z = — 1). Damped beatings at Sl± = 0±5 are possible 
due to the hybridization of fl « 8 (Fig. HJa) inset) . 

We illustrate these features in Fig. 2] for v — e. The 
beatings visibility is reduced with increasing T, due the 
onset of the pure dephasing processes. For an unpolar- 
ized state, p z = 0, (cr-(t)) shows a intermediate behav- 
ior between the ones at p z — ±1 since at resonance all 
eigenvalues contribute (Fig. EJd)). Damping is strongest 
for p z = — 1, weakest for p z — I and intermediate for 
p z — 0. In fact, for e > A, preparation in the pure state 
Pz = +1 makes the impurity close to its ground state and 
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FIG. 4: Ke[{<T-(t))/{a-(Q))] at resonance £ = v = 3A for 
initial states (a) p z = 1: slow decay with IZe [A1/4], (b) p z — 
— 1: fast decay with 72.e[A2/3]. Parameters: T — (blue), 
T = 0.5A (red), T = 0.9 A (green), K — 0.1. 



non Gaussian and back-action effects due to a coherent 
bistable impurity. These may represent a ultimate limi- 
tation for solid state qubits even when single shot mea- 
surement schemes are available. Our analysis by chang- 
ing temperature, strain e and coupling v, may provide 
valuable insights to realistic scenarios where a wide dis- 
tribution of the parameters has to be considered [loj ]. 
The employed SBE represents a general effective model 
for complex physical baths awaiting specific microscopic 
description, as those typical of solid state nanodevices. 



less damped, while it is close to the excited state when 
p z — — 1 with strongest damping. 

In the last part of this Communication we present an 
alternative perspective, considering the qubit as a mea- 
suring device for a mesoscopic system described by the 
SBE. Remarkably, the qubit acts as a detector despite 
the absence of direct qubit-SBE inelastic transitions [19J , 
In fact, the pure dephasing coupling amounts to a "dis- 
persive", quantum non-demolition measurement regime 
for the qubit. Detection is feasible due to the qubit 
back-action on the SBE. This point of view is illustrated 
in Fig. [5l where the length of the Bloch vector in the 
x - y plane, |(<r_(t))|, acts as a sensitive detector of the 
mesoscopic system ("impurity") preparation. At reso- 
nance, the unpolarized state, p z — 0, is identified by 
beatings, Fig. E]Ja). These almost disappear for pure 
states, p z = ±1, where oscillations at fi+ occur, Fig.OJb). 
Identification of the impurity preparation far from res- 
onace results instead from different oscillation amplitudes 
and/or decay rates, Fig. [5] (c) - (d). 

In conclusion, we have identified in time domain 
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FIG. 5: I <tr_ (*)> / <cr_ (0)> ] for e = 3A, K = 0.1. Panels (a) 
and (b): resonant impurity v — e. (a) p z = at T = (black) 
and T = 0.5 A (red); (b) T = for p z = 1 (orange), p z = -1 
(black) . Panels (c) and (d) : non resonant case v — A at T — 
(blue) and at T = 0.9A (red). In (c) p z = 0, in (d) p z = 1 
top, p z — — 1 bottom. Note the weak T-dependence. 



We acknowledge support from the EU-EuroSQIP (IST- 
3-015708-IP) and MIUR-PRIN2005 (2005022977). 



[1] 
[2] 
[3] 



[4 

[5. 

[6 

[7; 

[8 
[9 
[10 

[11 

[12. 

[13 

[i4; 

[15 
[16 

[17; 
[is: 

[19 



M.B. Weissman, Rev. Mod. Phys. 60, 537 (1988). 

A.B. Zorin et ai, Phys. Rev. B 53, 13682 (1996). 

T. Duty et al, Phys. Rev. B 69, 140503(R) (2004); M.G. 

Peters et ai, J. Appl. Phys. 86, 1523 (1999); J. Eroms et 

al, Appl. Phys. Lett. 89, 122516 (2006). 

R.W. Simmonds et al, Phys. Rev. Lett. 93, 077003 

(2004); K.B. Cooper et ai, ibid. 93, 180401 (2004). 

O. Astafiev et al, Phys. Rev. Lett. 93, 267007 (2004); 

0. Astafiev, et ai, Phys. Rev. Lett. 96, 137001 (2006). 
R. de Sousa et al, Phys. Rev. Lett. 95, 247006 (2005); 
L. Faoro, L.B. Ioffe, Phys. Rev. Lett. 96, 047001 (2006). 

1. Martin, Y. M. Galperin, Phys. Rev. B 73, 180201(R) 
(2006). 

A. Shnirman et al, Phys. Rev. Lett. 94, 127002 (2005). 
E. Paladino et al, Phys. Rev. Lett. 88, 228304 (2002). 
Y.M. Galperin, B.L. Altshuler, D.V. Shantsev, 
|cond-mat/0 3 12490 Y.M. Galperin et al, Phys. Rev. 
Lett. 96, 097009 (2006); J. Bergli, Y.M. Galperin, and 

B. L. Altshuler, Phys. Rev. B 74, 024509 (2006). 

G. Falci et al, Phys. Rev. Lett. 94, 167002 (2005). G. 

Ithier et ai, Phys. Rev. B, 72, 134519 (2005). 

Y. Nakamura, Yu.A. Pashkin, J.S. Tsai, Nature 398, 

786 (1999); D. Vion et al, Science 296, 886 (2002); 

J. Martinis et ai, Phys. Rev. Lett. 89, 117901 (2002); 

I. Chiorescu et al, Science 299, 1869 (2003). 

G. Falci et ai, Phys. Rev. A 70, R40101 (2004). 

U. Weiss Quantum Dissipative Systems 2nd Ed (World 

Scientific, Singapore 1999). 

E. Paladino, M. Sassetti, G. Falci, Chem. Phys. 296, 325 
(2004); E. Paladino et al, Chem. Phys. 322, 98 (2006). 
Numerical analysis of relaxation in H. Gassmann and C. 
Bruder, Phys. Rev. B 72, 035102 (2005). 
A. Grishin, I.V. Yurkevich, I.V. Lerner, Phys. Rev. B 72, 
060509(R) (2005). 

C. Cohen- Tannoudji, J. Dupont-Roc and G. Grynberg 
Atom-Photon Interactions, Wiley-Interscience (1993). 
R. Aguado and L. P. Kouwenhoven Phys. Rev. Lett. 84, 
1986 (2000). 



